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Following the idea of the density functional approach, we develop a generalized Bogoliubov theory 
of an interacting Bose gas confined in a one-dimensional harmonic trap, by using a local chemical 
potential - calculated with the Lieb-Liniger exact solution - as the exchange energy. At zero tem¬ 
perature, we use the theory to describe collective modes of a finite-particle system in all interaction 
regimes from the ideal gas limit, to the mean-field Thomas-Fermi regime, and to the strongly inter¬ 
acting Tonks-Girardeau regime. At finite temperature, we investigate the temperature dependence 
of collective modes in the weak-coupling regime by means of a Hartree-Fock-Bogoliubov theory with 
Popov approximation. By emphasizing the effects of finite particle number and nonzero temperature 
on collective mode frequencies, we make comparisons of our results with the recent experimental 
measurement [E. Haller et al, Science 325, 1224 (2009)] and some previous theoretical predic¬ 
tions. We show that the experimental data are still not fully explained within current theoretical 
framework. 

PACS numbers: 67.85.-d, 02.60.Cb 


I. INTRODUCTION 

Many-particle systems in one dimension (ID) with a 
short-range interparticle interaction plays an important 
role in understanding fascinating quantum many-body 
physics For instance, a ID interacting Bose gas 

is theoretically anticipated to experience different phases 
by changing interaction strength and temperature, and 
to exhibit a number of intriguing phenomena, such as 
effective fermionization and nontrivial quench dynam¬ 
ics P, Pi . Since the realization of Bose-Einstein conden¬ 
sation in 1995, more and more experiments have dealt 
with ID atomic bosons in a harmonic trap at ultra¬ 
cold temperature, which can help test theoretical pre¬ 
dictions and understand all these unusual phenomena. 
At present, many experiments have been conducted and 
controlled by means of Feshbach resonances P or di¬ 
rectly heating ID quantum degenerate Bose gases P, 
which measured and characterized variou^hysical quan- 
titieSjincluding momentum distribution P, pair correla¬ 
tion P and quenching rate 0, P- Motivated by these 
rapid experimental advances, there are numerous the¬ 
oretical studies based on sum-rule approach [91, y aria- 
tional method [ip , local density approximation [ll| and 
diffusion Monte Carlo simulations [ip- Experimental 
measured quantities, particularly momentum distribu¬ 
tion [ip and pair correlation [l3 - [lp , have been predicted 
and compared with experimental data. 

In the recent experiment P, the Feshbach resonance 
technique is used to tune the interatomic interaction of 
a ID Bose gas in a harmonic trap at extremely low 
temperature. The measured ratio of squared breath¬ 
ing mode frequency exhibits a reentrant be¬ 

haviour, from 4 in the non-interacting regime to 3 in 
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the weakly interacting regime, and then back to 4 in 
the Tonks-Girardeau regime [l3|. Most recently, this 
interesting reentrant behaviour was addressed by two 
theoretical works based on simulations at zero tempera¬ 
ture [nl, [Tpl . Choi and coworkers used a time-dependent 
modified nonlinear Schrodinger equation (m-NLSE) with 
a local chemical potential replacing the conventional non¬ 
linear term [HI- Gudyma and collaborators combined a 
sum-rule approach in the mean-field regime and the local 
density approximation in the Tonks-Girardeau regime to 
describe the breathing mode [Ip . Diffusion Monte Carlo 
simulations for few particle numbers were performed, in 
order to obtain density profiles as inputs to the sum-rule 
approach. Comparing all results together, we find that 
there are still deviations between experimental data and 
theoretical predictions, particularly in the deep weakly- 
interacting regime. 

By considering finite particle number (i.e., varying in 
a range of 8 ~ 25) and un-avoidable nonzero tempera¬ 
ture in the realistic experiment P, we therefore would 
like to address in a more systematic way their effects on 
collective mode frequencies of a ID trapped Bose gas. 
We note that, at large number of particles, the temper¬ 
ature dependence of mode frequencies has been recently 
investigated by using a hydrodynamic theory la¬ 
in this paper, we investigate low-lying collective mode 
frequencies of a harmonically trapped Bose gas in ID 
at a wide range of effective interaction parameter 7^^^, 
covering all interaction regimes from the non-interacting 
regime, mean-field regime to Tonks-Girardeau regime, by 
developing a generalized Bogoliubov theory at T = 0 and 
a Hartree-Fock-Bogoliubov th eory with Popov approxi¬ 
mation at finite temperature [iSj, Il9| . The former the¬ 
ory concentrates on the zero temperature case, where we 
follow the idea of the density-functional approach, use 
a generalized Gross-Pitaevskii equation and take the lo¬ 
cal chemical potential for uniform density obtained from 
the Lieb-Liniger model as the exchange energy (i.e., the 
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nonliner term) [^ . In the latter, the finite temperature 
effect is taken into account through the self-consistent 
Hartree-Fock-Bogoliubov equations in the weakly inter¬ 
acting regime. This is the only theory that we know so 
far to address both effects of nonzero temperature and 
finite particle number. 

This paper is organized as follows. We provide details 
of the two theoretical methods in Sec. im In Sec. imi 
we present the ground state of ID Bose gases (see Figs.[T] 
and[ 2 ]) and discuss the effects of particle number and tem¬ 
perature on excitation frequencies (see Figs. [310] and [5]). 
We then compare our results of the breathing mode fre¬ 
quency with the experimental data Q and the previous 
theoretical predictions [ll| (see Figs. ElandjT]). Mode fre¬ 
quencies of two higher collective modes are also discussed 
in Fig. El Summary and outlook are given in Sec. IlYl 


II. THEORETICAL FRAMEWORK 


A. A density functional method with Lieb-Liniger 
integrals at zero temperature 


We start with a ID atomic Bose gas with a repulsive 
zero-range potential. The system of N particles can be 
described with the Lieb-Liniger Hamiltonian [^ . 


7 y = -^V 

9.m 


2m ^ dxj 

i—1 ^ 


N 
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where m is the mass of atom and is the ID interaction 
strength of form 


9id 


2h^ 


( 2 ) 


characterizing the interaction between bosons, is the 
ID scattering length calculated by plj l 


2a, 


1 -C^ 

a_L 


( 3 ) 


where ajj, is the three-dimensional (3D) scattering 
length, a_L = [?i/(tow_l)]^/^ is a two-dimensional (2D) 
harmonic oscillator characteristic length (see below) and 
the constant C ~ 1.4603. Experimentally, systems of ID 
Bose gases are usually trapped in a cylindrically sym¬ 
metric potential with strong transverse confinement and 
weak longitudinal confinement (i.e., with trapping fre¬ 
quencies a;_L ^ Wx) • Transverse excitations are not taken 
into account if the transverse vibrational energy Sajj_ is 
much larger than the chemical potential or the energy 
scale of the thermal cloud (fiw_L » /r, k^T) How¬ 

ever, the scattering of two atoms in the lowest transverse 
mode could strongly be affected by high transverse ex¬ 
citations, when the 3D scattering length 031 , is close to 
the 2D harmonic oscillator length a_L. This leads to a 
confinement-induced Feshbach resonance, as can be seen 
in Eq. (|3]). 


1. Lieb-Liniger model 


Lieb and Liniger investigated this system at T = 0 and 
solved it exactly in 1960s. They dealt with the Hamil¬ 
tonian (ED and the associated Schrbdinger equation by 
means of Bethe ansatz (2^, and obtained exactly the 
ground state as well as low-lying excited states in a uni¬ 
form gas with a constant density n = NjV [^ . 

In the Lieb-Liniger model, they defined a dimensionless 
interaction parameter 7 (n) as a function of , which is 


7(n) 


"T-glP 


2 
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By solving exactly the Hamiltonian, they got a group 
of integral equations (i.e., Lieb-Liniger integrals), includ¬ 
ing the normalization condition and the equation of the 
ground state energy. The ground state energy has the 
form 


Eo = 


Nh^n^ 

2m 


e(7), 


( 5 ) 


where 0 ( 7 ) is a dimensionless function of 7 , which can be 
obtained by solving the following Lieb-Liniger integrals. 


9{x) = 


2 \ 


2tt 27r {x — x')"^ 

^ = J 9 ix)dx, 

e(7) = {jf J ^x^9{x)dx. 


g{x')dx', 


( 6 ) 


In the integrals, g{x) is the distribution function of the 
quasi-momentum, which is represented by the variable 
X after a rescaling Thus, the variables x,x' are all 
bounded in the range [— 1 , 1 ], as shown in the upper and 
lower limits of the integrals. Variable A is proportional to 
7 since the remaining integral part in A is fixed. Hence, 
A is also proportional to the interaction strength for 
a fixed uniform density n. 

By numerically solving the integrals, it can be shown 
that 6 ( 7 ) is a monotonically increasing function of 7 [ 2 ^. 
In the limit of 7 = 0, e(0) = 0 and the ground state 
energy is Eq — 0 , corresponding to the case that all free 
bosons occupy the zero-momentum state. When 7 is suf¬ 
ficiently large, the asymptotic value of e(j) is 7 r^/ 3 , which 
exactly coincides the predicted value for impenetrable 
bosons in ID by Girardeau El- All the information of 
the ground state of ID Bose gas can then be calculated 
from 6 ( 7 ). The energy per particle at the ground state is 
e(n) = Eq/N = h‘^n'^e{'y)/{2m), and the corresponding 
chemical potential /i(n) is given by 


H{n) = 


d{ne{n)) 


dn 


2m 


9il)- 


( 7 ) 


It is straightforward to see that the dimensionless chem¬ 
ical potential ^( 7 ) can be calculated by 


^( 7 ) = 3e(7) - 7e'(7). 


( 8 ) 
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In Ref. [HI, Choi and coworkers proposed the following 
analytic expression for the chemical potential at large 7 , 


fi{n) 


7 ^(n )(2 + 87 ( 71 )) 

2 m _ 3(2 + 7 ( 71 ))^ _ ’ 


(9) 


near the Tonks-Girardeau limit. Compared with the nu¬ 
merical results of Lieb-Liniger integrals, this expression 
turns out to be very accurate, as long as 7 is larger than 
10 . 


2. Generalized Gross-Pitaevskii theory 

The Gross-Pitaevskii theory provides a good descrip¬ 
tion of weakly-interacting atomic Bose-Einstein conden¬ 
sates at T = 0. However, it is known that the theory 
fails in low dimensions {d < 2) and therefore needs ap¬ 
propriate modifications. In Ref. [ 2 ^, Kolomeisky and 
coworkers suggested an interesting modification in the 
Tonks-Girardeau limit for ID trapped Bose gases. Here, 
we generalize their idea to all interaction strengths, fol¬ 
lowing the procedure of using an improved exchange en¬ 
ergy in the density-function approach. 

Recall that the standard Gross-Pitaevskii equation 
with a condensate order parameter (wavefunction) 
'I'(a;,t) is given by 

= [^-0 +giD'n{x,t)]'it{x,t) ( 10 ) 

with the single particle Hamiltonian 

where 14 xt(a;) = harmonic trapping po¬ 

tential and 77(0;, t) = I (x, t)p is the particle density. The 
Hartree term g-^^n{x,t) in Eq. (fTUl) is only applicable in 
the weak coupling regime. 

Away from the weak coupling regime {na^-^ ^ 1 or 
7 <C 1), we may use local density approximation (EDA) 
to determine the ground state of a trapped system. The 
LDA amounts to setting, 

/i — /Xloc(77(x)) -t- V^xt(^)) (^^) 


Modifying the interaction term to the local chem¬ 
ical potential ii\oc(n), the generalized Gross-Pitaevskii 
equation reads, 

^ (13) 

which has the same form as the m-NLSE equation 
adopted by Choi and coworkers m- 

The idea of directly modifying the exchange-energy¬ 
like term at zero temperature is supported by the follow¬ 
ing derivation of hydrodynamic equation. That is, we 
rewrite the order parameter as 

«'(a;,t) = 775 ( 1 , (14) 


where n{x, t) is now interpreted as the time-dependent 
superfluid density of the system and 0{x, f) is the associ¬ 
ated phase. We therefore introduce a superfluid velocity 
field 'v(x,t) = ^ and rewrite the generalized 

Gross-Pitaevskii equation in terms of the superfluid den¬ 
sity n{x,f) and superfluid velocity v(x,<): 


9V d , , , rr / N 1 

m— + —[pioc[n) + Ve^t[x) + -mv ) = 0 , 

dn d , . ^ 

_ + _(„v) = 0. 


(15) 


Here we have neglected a quantum pressure term, which 
is small in the long wave-length limit. Thus, it is clear 
that our generalized Gross-Pitaevskii equation is iden¬ 
tical to the standard ID time-dependent hydrodynamic 
equations in the long wave-length limit. For a finite parti¬ 
cle number system, it is preferable to use the generalized 
Gross-Pitaevskii equation to understand the dynamics of 
interacting ID Bose gases. The effect of finite particle 
number is taken into account by the quantum pressure 
term that is neglected in the hydrodynamic equations. 


3. Generalized Bogoliubov theory 

We now consider excited states of ID Bose gases, which 
can be treated as small oscillations around the superfluid 
density at the ground state (i.e., foix)) with frequencies 
coi- Their wavefunctions, Ui and Uj, are given by [ 2 ^, 


with fj, being the global chemical potential. Once we 
know the local chemical piociji) for a uniform Bose gas 
from the Lieb-Liniger model, we can then determined 
the density profile n{x) by inversely solving Eq. (fT^ .The 
important observation made by Kolomeisky and cowork¬ 
ers is that in the Tonks-Girardeau limit, one may sim¬ 
ply obtain a modified Gross-Pitaevskii equation by us¬ 
ing g, ]„r.( n) = /2m to replace the Hartree term 

5id^ [23[. Motivated by this work, one can use isioc{n) 
calculated numerically in the Lieb-Liniger model, to ob¬ 
tain the corresponding result in the intermediate regime 
between the mean-field limit and the Tonks-Girardeau 
limit. 


^(x,t) = \Mx) + u^{x)e-^‘^^*/’^ + v*{x)e^"^*/^ 


The corresponding density 77 


77 = |(^o|" 




= |'I'(x,t)p is 
-F -f H.c. 


(16) 

(17) 


where we keep only the linear terms in the complex func¬ 
tions Ui and Vi. Accordingly, the local chemical potential 
Mioc( 77 ) can be written as 


Mloc(77) 


Woe ( 770 ) + 


dp\QC 

dn 


+ H.c. 
(18) 
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in the Taylor expansion with no{x) = \(j)o(x)\'^, where 
higher orders are neglected. Taking Eqs. (HID and m 
back into Eq. (fT^ and sorting out the terms in e 

^ obtains respectively the 
static generalized Gross-Pitaevskii equation, 

^(/)o(a;) = [Ho + ^J.loc{no)] 0o(a;), (19) 


as well as the coupled Bogoliubov equations 


CUi{x) + AiVi{x) = + fiLUiUi{x), 
Cvi{x) + Mui{x) = — hwiVi{x), 


( 20 ) 


where we have defined the operators 

dnioc{n) 


C — 'Ho ~ fJ- H fJ-iocino) + 


dn 


no, (21) 


j n—UQ 


and 


M = 


^Mioc(n) 

dn 


no- 


( 22 ) 


The Bogoliubov wave-functions Ui{x) and Vi{x) satisfy 
the normalization condition, 


/ OO 

dx{u*{x)uj{x) — v*(x)vj(x)) = 6ij. (23) 

-CX) 

I 

= J dx^'^{x,t){Ho - fJ,)^ix,t)+ ^ JJ 


where '^{x,t) is the bosonic field operator. By taking a 
(5-function interaction U{x — x') = g-^^5(x — x') between 
bosons, the Heisenberg equation of motion is given by 


ih 


5^(x,t) 

Wt 


(Ho - M)4'(x,t) 


We separate the operator y(x, t) into the condensate part 
and non-condensate part |25l l30l - l32j : 


^{x,t) = ^{x,t) + fi{x,t), (26) 

which leads to ('I') = $ and (?)) = 0. Thus, we obtain 
(^I'h) = -I- ( 77 I 77 ) and (^T) = $^ -|- {fjfj). 

Within the mean-field approximation, we rewrite the 
term as [l^ 

(^t^,^) = _|_ 2(f)t^))$ -I- (r)r))$*. (27) 


The above formalism (with fj-iocin) = g^j^n) was first 
introduced by Pitaevskii, in order to investigate excita¬ 
tions of a vortex line in a uniform Bose gas. One can 
get the same result if one diagonalizes the Hamiltonian 
with Bogoliubov transformation 2^. In 1996, Burnett 
and colleagues used the coupled equations (|2nil to study 
the properties of excited states in 3D weakly-interacting 
trapped Bose gases [ 2 ^. 

After numerically solving the static Gross-Pitaevskii 
equation (HU) and the coupled equations (EOl) , one can 
obtain directly the profile of the ground state, as well as 
excited states with energies = hui and hence the fre¬ 
quency Wni of the breathing or compressional monopole 
mode. 


B. Hartree-Fock-Bogoliubov method with Popov 
approximation at finite temperature 

In quantum many-body systems, the Hartree-Fock ap¬ 
proximation (or self-consistent approximation) provides 
a good description of the ground state wavefunction and 
energy at weak couplings. It works for both bosons and 
fermions. It can also be generalized to finite temper¬ 
ature to solve thermodynamics [ 2 ^. For a 3D weakly 
interacting Bose gas, Shi and Griffin have shown that 
the Hartree-Fock theory is useful at finite temperature 
except a small critical area near the transition tempera¬ 
ture [H,!!!. 

In the case of a ID trapped Bose gas, the grand- 
canonical Hamiltonian takes the form 


dxdx''^^ (x,t)'^^ {x',t)U(x — x')'^{x' ,t)'^{x,t), (24) 


By taking an average on both sides of Eq. (EH) and mak¬ 
ing use of Eq. (EZl), we then find that 


94? 

ih— = {Ho - + 5id [{nc + 2nt)$ -b na$*], (28) 


where Uc = |4)(x,t)p is the time-dependent condensate 
density, nt = is the non-condensate thermal density 
and Ua = {gg) is the anomalous thermal density. 

The equation of motion for the non-condensate opera¬ 
tor g can be similarly obtained: 


dfj 

ih— = {Ho - g)g + SiD [2(nc + nt)g + (ric + «a)?7^] ■ 

(29) 

To solve this equation, we then use the Bogoliubov trans¬ 
formation to expand the non-condensate operator g{x, t) 
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and its conjugate as 

i 

f]\x,t) = e*^*/'*^ \vi{x)e~'‘‘^''^oti + u*{x)e'‘‘^'*a\ , 

(30) 

i.e., the non-condensate operator fj{x,t) is rewritten in 
a quasi-particle basis. The operators d^ and a are the 
creation and annihilation operators for quasi-particles 
respectively and, Ui and vt are the corresponding am¬ 
plitudes. Bogoliubov quasi-particles are assumed to be 
non-interacting, and their operators d^ and d satisfy the 
bosonic commutation relations: 


a\ , ttj 


— 6ij, 


d( d( 


— [d^, cyj ] — 0 . 


(31) 


Taking the Popov approximation, neglecting the 
anomalous term with Uq = (r)^) we finally obtain 
the time-dependent Hartree-Fock-Bogoliubov equations 
(HFBP) [la 112) [13 j which consists of two parts: 

(A) a modihed Gross-Pitaevskii equation 

5$ 

ih—= [no - ^J, +g^^{nc + 2nt)]^, (32) 

which leads to 


g(j){x) = [Ho + g^oinoix) + 2nt{x))] (j)[x,t)- (33) 

Here $(a:,t) = and fl = + fj,. We have 

introduced a small chemical potential difference A/i to 
allow a finite particle number of condensate (see Eq. dSZD 
below). The static condensate density is no = \4‘{x)\^ 
and the thermal density is nt(x) = 

(B) coupled Bogoliubov equations 


■ £ 

M ' 

Ui{x) 

= hwi 

+Ui{x) 

M 

C 

_ Vi{x) _ 

_ -Vi{x) _ 


which are obtained by substituting Eq. (15(11) back into the 
equation of motion for the non-condensate operator t). In 
the coupled equations, we have defined 

C = no-ij- + 2g^^{no{x)+nt{x)), (35) 


and 


M=g^j,no{x). (36) 


In the modified Gross-Pitaevskii equation, there are 
condensate and non-condensate densities in the exchange 
term, both of which have to be solved self-consistently. 
It should be noted that, for the condensate density, be¬ 
cause our bosonic system has a finite particle number, 
we have to introduce a small chemical potential differ¬ 
ence to account for the finite condensate particle number 
Nq = f dxnoix), i.e.. 


No = 


I 

g/3A/i _ I ’ 


(37) 


with the inverse temperature /3 = l/k^T. The quasi¬ 
particle occupation number is affected by the finite con¬ 
dense number as well. For the i-th quasi-particle occu¬ 
pation number Ni = (djdi), we have, 

pP{hu)i+Afi) _ 1 

1 (38) 

= (i + ^)e/3fi-i-I- 
Therefore, the thermal density is given by 
nt{x) = {f|^) 

= +«f(^)] (1 ^ - 1 + 

= E 

i i 

(39) 

The chemical potential of the system, /i, is to be deter¬ 
mined by the number equation for the total number of 
atoms, 

N = No+Nt, (40) 

where Nt = f dxnt{x). 


III. RESULTS AND DISCUSSIONS 

We are now ready to perform numerical calculations of 
the above mentioned approaches and compare our results 
with the experiment data Q and the previous theoretical 
predictions m- In the experiment, a 2D optical lattice is 
used to trap about (I ~ 4) x 10“^ Gs atoms in (3 ^ 6) x 10^ 
ID tubes with 8 ^ 25 atoms in the center tube. Ghoi 
et al. dealt with the case of a particle number N = 
25 and introduced an effective dimensionless interaction 
parameter , which is defined as 

2 Q,T,7r 

--qlq- 1 = %7' 

^tg(^)I^idI '\/ 2 N 

Here, n.p( 3 ( 0 ) = 2NmuJho/fi/ir is the peak density in 

the Tonks-Girardeau regime at the tube center [^ . 

In our calculations, harmonic oscillator units are used 
with h = Who = m = I and fee = I. Length and 
energy are written in the units of harmonic oscillator 
length Oho = [^/(jTiWho)]^^^ and harmonic oscillator en¬ 
ergy fwjho = respectively. 

We fix the particle number at A = 8,17, 25, and vary 
the effective interaction parameter . The ground state 

and collective modes of the ID harmonically trapped 
Bose gas at T = 0 are obtained, by numerically solv¬ 
ing the generalized Gross-Pitaevskii equation m, and 
the generalized Bogoliubov equations O, where the lo¬ 
cal chemical potential is obtained with Lieb-Liniger inte¬ 
grals [ 13 . Then we introduce the non-condensate ther¬ 
mal density nt by means of the HFBP theory, and com¬ 
pare the ground states and collective behaviours at finite 
T with the result at T = 0. 
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FIG. 1. (Color online) Density profiles in the mean-field 
regime (a) and in the Tonks-Girardeau regime (b), where the 
effective interaction parameter = 10”^ and 10 have been 
used, respectively. The results are calculated with N = 8 
(black solid lines), = 17 (red dashed lines) and N = 25 
(blue dotted lines). The density n(x) and position x are taken 
in units of harmonic oscillator length = ■s/h/(rnuJho) and 
^ho’ respectively. 


FIG. 2. (Color online) Density profile: at T = 0 predicted by 
the generalized Gross-Pitaevskii theory (black circles); at T = 
0 (red dashed lines) and hnite temperatures (other colorful 
lines) calculated by the HFBP theory. The results are shown 
at the effective interaction parameter = 10“^ in the mean- 
field regime. The particle number N is hxed at 25. is the 
critical temperature for a ID ideal Bose gas, which can be 
estimated as ksT^ = fkJhoX/ ln(2N) (^. 1^ . 


A. Density profile 

In this subsection, we study the density profile of the 
ground states. In particular, at zero temperature we 
compare the results obtained by the generalized Gross- 
Pitaevskii theory with those predicted by the HFBP the¬ 
ory. 


presented by black circle for comparison. It agrees well 
with the prediction of the HFBP theory (i.e., the red 
dashed line). With increasing temperature, the conden¬ 
sate fraction decreases, as well as the condensate occu¬ 
pation number, leading to the decreasing of the height of 
density profiles. 


B. Breathing mode 


1. T = 0 case/ the generalized Gross-Pitaevskii theory 

The ground states for different particle number N = 
8,17, 25 at T = 0 are shown in Fig.[Tl obtained by solving 
the generalized Gross-Pitaevskii theory. In the figure, we 
focus on the mean-held regime (with = 10“^, the left 
panel) and the Tonks-Girardeau regime (with y^^^ = 10, 
the right panel). At y^^j = 10“^, the density prohle is 
roughly a Gaussian curve, while at y^jj = 10 it tends to be 
a semicircle. In both regimes, for a given y^^j, the height 
of density prohles is enhanced with increasing particle 
number. 


2. Finite T at weak couplings: the HFBP theory 

The ID HFBP theory is valid only in the weak cou¬ 
pling regime. The corresponding density prohles at dif¬ 
ferent temperatures are shown in Fig. [2] For T = 0, 
the density prohle predicted by the generalized Gross- 
Pitaevskii theory (indicated as gGP in the hgure) is also 


In this subsection, collective modes, especially the 
breathing mode, are investigated with the generalized 
Bogoliubov theory and the HFBP theory. The dipole 
mode frequency should precisely be the trapping fre¬ 
quency Who, according to the Kohn theorem. We recover 
the result with a relative error about 0.1% ^ 0.2% with 
respect to Who- 


1. Particle number effect in two limits 

The experiment on the breathing mode frequency was 
conducted for particle numbers in the range of 8 ^ 25. 
There is a deviation between the experimental results Q 
and one of previous numerical simulations at A = 25 in 
the weak coupling regime m- In order to check whether 
this is due to the effect of different particle number, we 
calculate the breathing mode frequency with different N. 

In the deep weak coupling limit, it is convenient to 
define a Hartree parameter A = |/(fVaho) ^ 1- The 
sum-rule approach predicts that in the limit of A ^ 1, 
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FIG. 3. (Color online) The ratio of the squared breathing 
mode frequency otm/iUho as a function of the particle number 
N. Our numerical calculation is shown by the black solid line, 
and the analytic result Eq. (H21) is shown by the red dashed 
line. All results are near the non-interacting limit with an 
interaction strength = 10“® in the trap units. 


FIG. 4. (Color online) The ratio of the squared breath¬ 
ing mode frequency as a function of C,^y/N (black 

squares), compared with the analytic prediction (Eq. (1431) ') 
that is shown by the red dashed line. The inset shows the 
same ratio as a function of the particle number N. Here, we 
take an interaction strength = 10^. 


the correction of a finite particle number on the squared 
breathing mode frequency is [T^ 

;|l:^4(l-c,A-i) 

(42) 

= 4_^iV,A^oo 


where = l/\f^ for all > 2. This analytic predic¬ 
tion indicates that the squared breathing mode frequency 
ratio is proportional to the particle number N 

at a fixed . 

In Fig. [21 we show our numerical results (solid line) 
at A = |a^jj |/(A^aho) 1, with the effective interaction 

parameter varying from 4.4 x 10“^ to 1.6 x 10“^ at 
a constant = 10“^. The analytic results Eq. (l42ll 
are also shown by the dashed line. There is a very good 
agreement, within a relative error 0.1%. Presumably, 
this small discrepancy is due to the over-estimation of the 
mode frequency in the sum-rule approach, which predicts 
only an upper bound for the mode frequency 371. 

In the opposite limit of strong couplings, which is char¬ 
acterized by the parameter A = /alo <C 1, the cor¬ 
rection on the breathing mode frequency due to a finite 
N is also known: [M 


^=.4(1-C,v/A) 

= 4-—C'n\/]V,A^0 

^ID 

where is given for all > 2 by 

-\)T{N+ \) 

^ ~ TT^ r{N)r{N + 2) 

X3 ^2 0, 1 - N, -TV; I - iV, i - iV; 1^ . 


(43) 


(44) 


The Af-dependence of C,^ is very weak and its value varies 
from Cj Ri 0.282 to fv 0.306, and to ~ 0.305 in 
our case. It is clear from Eq. (H21) that for a fixed the 
ratio of the squared breathing mode frequency 
has a linear dependence on the combined variable vV. 

In Eig. m we verify this linear behaviour by taking a 
fixed g^j^ = 10^ in the Tonks-Girardeau regime, for which 
the effective interaction parameter varies from 44 to 
157. 


8. Finite temperature effect in the weak coupling limit 


Here we consider the finite temperature effect in the 
weakly interacting limit, by calculating collective mode 
frequencies using the weak-coupling HEBP theory. It is 
known that at sufficiently low temperature, this effect is 
small since the thermal fraction of the system is negligi¬ 
ble.For example, Debbie Jin’s group has shown that the 
measured collective oscillating frequencies of a 3D Bose 
gas at temperature T < 0.48T° have a good agreement 
with the theoretical predictions at T = 0 [s^. To em¬ 
phasize the effect of finite temperature on the collective 
mode, we consider here T > 0.48T°. 

In Fig. 121 we compare the ratios of the squared breath¬ 
ing mode frequency at T = 0 and T = 0.8r°. The mode 
frequency becomes larger at finite temperature. At suffi¬ 
ciently large temperature, actually we anticipate that the 
ratio approaches the ideal gas limit, i.e., = 4. We 

also compare the zero temperature ratios, predicted by 
the generalized Bogoliubov theory (symbols, indicated 
as gB in the figure) and the HFBP theory (dashed line). 
There is a good agreement, as anticipated. 
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FIG. 5. (Color online) The ratio of the squared breathing 
mode frequency at = 25 in the weakly-interacting 

regime. We have shown the results at T = 0 predicted by 
the generalized Bogoliubov theory (black squares) and by the 
HFBP theory (red dashed line) and the results at T = 0.8T° 
given by the HFBP theory (blue dotted line). 



3. Comparisons with the experiment and previous theory 


We now address the breathing mode frequency in all 
interaction regimes, emphasizing its dependence on fi¬ 
nite particle number and nonzero temperature. We vary 
the effective interaction parameter 7 ,.^^ from 2.2 x 10 “^ 
to 2.9 X 10^, and thus cover all the regimes from the 
non-interacting limit, the mean-field regime to the Tonks- 
Girardeau limit. The results are presented as a function 
of the interaction parameter shown in the linear (Fig. [ 6 ]) 
and logarithmic scales (Fig. [7]), in comparison with the 
experimental data Q and a previous theoretical predic¬ 
tion (TTj |. 

In general, in these figures the squared frequency ratio 
^m/^ho of breathing mode decreases from 4 in the 
non-interacting limit to 3 in the weakly-interacting mean- 
field regime, and then increase back to 4 again in the 
strongly-interacting Tonks-Girardeau regime. 

In greater detail, the previous theoretical work (see the 
results indicated as m-NLSE in the figures) considered a 
particle number N = 25 ll|. Here we have performed 
numerical calculations with the same number of parti¬ 
cles. We have also considered other two sets of particle 
number, N = 8 and N = 17, since in the experiment the 
range of the particle number iV is 8 ^ 25 Q . For the case 
with a particle number N = 25, our results agree very 
well with the m-NLSE predictions. The good agreement 
is easy to understand, as both theories start from the 
same generalized Gross-Pitaevskii equation. The differ¬ 
ent numerical treatments, i.e., the time-dependent sim¬ 
ulations in Ref. m and our solution of the general¬ 
ized Bogoliubov equations, only lead to a negligible dif¬ 
ference. By further comparing both theoretical predic¬ 
tions at = 25 with the experimental data, we find a 
good agreement in the mean-field and Tonks-Girardeau 
regimes, where the breathing mode frequency essentially 
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FIG. 6 . (Golor online) The ratio of the squared breathing 
mode frequency Wm/Who as a function of the effective inter¬ 
action parameter covers all interaction regimes and 

varies from 2.2 x 10“® to 2.9 x 10^. We consider three par¬ 
ticle numbers: N = 8 (black solid line), N = 17 (red dashed 
line) and N = 25 (blue dotted line). We have compared our 
result with a previous theoretical prediction obtained by us¬ 
ing time-dependent modified nonlinear Schrodinger equation 
(m-NLSE) (yellow dot-dashed line) [lj| and the experimental 
data (green squares with error bars) [^. 
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- 2-1012 


FIG. 7. (Golor online) The ratio of the squared breathing 
mode frequency 8 -s a function of logj^g 7 ^^^. The plot is 

the same as Fig.[ 6 l but is shown here as a function of the inter¬ 
action parameter in a logarithmic scale, in order to emphasize 
the particle number dependence in the non-interacting limit. 


does not depend on the particle number. However, near 
the non-interacting limit, the discrepancy between ex¬ 
periment and theory becomes evident: the experimental 
data lie systematically below the theory curves. In this 
limit, the particle number dependence of the breathing 
mode frequency is significant. 

The particle number dependence is particularly clearly 
seen in Fig. [T] The decreasing of the particle number N 
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FIG. 8. (Color online) The frequency of higher order compres- 
sional modes, wsrd/ti^ho (a) and ct; 4 th/ciJho (b), as a function of 
logio7effi 8-t different particle numbers: N = 8 (black solid 
line), = 17 (red dashed line) and N = 25 (blue dotted 
line). 


from 25 to 8 increases the ratio of the squared breathing 
mode frequency. Thus, taking into account the possibil¬ 
ity of a smaller particle number (i.e., N < 25) in the real 
experiment will even enlarge the discrepancy between ex¬ 
periment and theory. On the other hand, this discrep¬ 
ancy cannot be resolved as a finite temperature effect, 
as in the previous subsection we have already examined 
that a nonzero temperature generally leads to a larger 
mode frequency. 


C. Higher order collective modes 

One of the advantages of our generalized Bogoliubov 
theory is that we can directly obtain higher order col¬ 
lective mode frequencies from numerical calculations. In 
Fig. [51 we present the mode frequency of the 3-rd (lower 
panel) and 4-th modes (upper panel) for all interaction 
regimes at T = 0, as a function of the effective inter¬ 
action parameter log^g . Similarly, the frequencies of 
higher modes exhibit the same reentrant behavior as the 
breathing mode frequency. 

Three analytic results can be used to understand the 
reentrant behavior. In the non-interacting limit, the 
mode frequency of the n-th mode is simply nwho- In 
the mean-field regime with sufficiently large number of 
particles, the collective mode frequency can be analyt¬ 
ically determined from a hydrodynamic theory, which 
predicts a;„ = ^n{n + l)/2a;ho [1, [13 ■ Therefore, we 
have Wg ^ 2.45a;ho and ~ 3.16who if iV —?> oo. With 
increasing number of particles, the minimum mode fre¬ 


quencies shown in Fig. [8] seem to approach these limiting 
values. Finally, in the Tonks-Girardeau limit, the mode 
frequency of the n-th mode again approach nwho, due to 
the effective fermionization of the system 

It is interesting to note that for high-lying collective 
modes, the effect of a finite particle number also becomes 
significant in the Tonks-Girardeau regime as well as in 
the mean-field regime. This is particularly evident for 
the 4-th mode, as shown in the upper panel of Fig. |8l 


IV. SUMMARY AND OUTLOOK 


In conclusions, we have studied collective modes of a 
ID harmonically trapped Bose gas, by developing two 
numerical approaches: (a) a generalized Bogoliubov the¬ 
ory at zero temperature following the idea of the density- 
functional approach, where we have taken the local chem¬ 
ical potential calculated from the Lieb-Liniger model as 
the exchange energy; (b) a conventional Hartree-Fock- 
Bogoliubov theory with Popov approximation at finite 
temperature, where the thermal density is included to 
make the theory self-consistent. By using these two ap¬ 
proaches and by emphasizing the effect of finite particle 
number and nonzero temperature, we have presented a 
systematic investigation of the breathing mode frequency 
in all interaction regimes and have explained the reen¬ 
trant behaviour of the mode frequency, which varies from 
4 in the non-interacting limit, to 3 in the mean-field 
regime and then back to 4 again in the Tonks-Girardeau 
limit. The frequency of higher order collective modes 
exhibits a similar reentrant behaviour. 

We have compared our result with the recent experi¬ 
mental measurement and a previous theoretical predic¬ 
tion. While our result agrees well with the previous theo¬ 
retical prediction, we have found that both theories can¬ 
not explain the measured mode frequency in the non¬ 
interacting limit. The discrepancy between experiment 
and theory becomes even larger when we take a small 
number of particles or a non-zero temperature. There¬ 
fore, we believe more theoretical investigations should be 
committed in the future in order to fully solve the dis¬ 
crepancy. Those works could focus on the issues such 
as the inter-tube tunneling and the confinement-induced 
three-body interparticle interaction. 
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